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ABSTRACT 

Searching for the signal of primordial gravitational waves in the B-modes (BB) power spectrum is one 
of the key scientific aims of the cosmic microwave background (CMB) polarization experiments. However, 
this could be easily contaminated by several foreground issues, such as the interstellar dust grains and the 
galactic cyclotron electrons. In this paper we study another mechanism, the cosmic birefringence, which can 
be introduced by a CPT-violating interaction between CMB photons and an external scalar field. Such kind of 
interaction could give rise to the rotation of the linear polarization state of CMB photons, and consequently 
induce the CMB BB power spectrum, which could mimic the signal of primordial gravitational waves at large 
scales. With the recently released polarization data of BICEP2 and the joint analysis data of BICEP2/Keck 
Array and Planck, we perform a global fitting analysis on constraining the tensor-to-scalar ratio r by considering 
the polarization rotation angle [ct(n)] which can be separated into a background isotropic part [a] and a small 
anisotropic part [Aa(h)]. Since the data of BICEP2 and Keck Array experiments have already been corrected 
by using the “self-calibration” method, here we mainly focus on the effects from the anisotropies of CMB 
polarization rotation angle. We find that including Aa(n) in the analysis could slightly weaken the constraints 
on the tensor-to-scalar ratio r, when using current CMB polarization measurements. We also simulate the mock 
CMB data with the BICEP3-like sensitivity. Very interestingly, we find that if the effects of the anisotropic 
polarization rotation angle can not be taken into account properly in the analysis, the constraints on r will be 
dramatically biased. This implies that we need to break the degeneracy between the anisotropies of the CMB 
polarization rotation angle and the CMB primordial tensor perturbations, in order to measure the signal of 
primordial gravitational waves accurately. 

Subject headings: cosmic microwave background - cosmological parameters - cosmology: theory 


1. INTRODUCTION 

Observations of the Cosmic Microwave Background 
(CMB) with high precisions are of great importance in mod¬ 
ern cosmology, since it is the most powerful way to under¬ 
stand the origin and evolution of our Universe. It is commonly 
believed that currently observed CMB anisotropies and large 
scale structure are originated from the quantum fluctuations 
at very early times preceding the hot expansion. Theoreti¬ 
cal hypothesis about this very early universe includes infla¬ 
tion (Guth 1981; Albrecht & Steinhardt 1982; Linde 1982), 
ekpyrotic/cyclic universe (Khoury et al. 2001) and so on. Ac¬ 
cording to these theories the Universe was smoothed and flat¬ 
tened at the early epoch, and at the same time the quantum 
fluctuations from vacuum extended outside the horizon and 
formed the primordial perturbations, i.e., the initial conditions 
of cosmic structure formation. 

The primordial perturbations can be decomposed into three 
types: scalar, vector and tensor perturbations. The scalar per¬ 
turbations cause the anisotropies of CMB and the struc¬ 
ture formation of large scale structure, while the vector 
perturbations, correspond to vortices decline rapidly dur¬ 
ing the expanding universe and can be neglected at later 
time (see e.g. Mukhanov (2005)). More interestingly, the 
tensor perturbations, i.e., the primordial gravitational waves 
(Starobinsky 1979; Rubakov, Sazhin & Veryaskin 1982; 
Fabbri & Pollock 1983; Abbott & Wise 1984) leave the 
unique imprint on the large scales power spectrum of CMB B- 
modes polarization (Kamionkowski, Kosowsky & Stebbins 


1997). Therefore, probes of CMB primordial B-modes po¬ 
larization are always treated as crucial ways to test the early 
universe models. Searching for the primordial gravitational 
waves is one of the key scientific goals of the CMB po¬ 
larization experiments, such as the BICEP3 (Ahmed et al. 
2014), the EBEx (Reichborn-Kjennerud et al. 2010), the SPT 
polarization (Hanson et al. 2013) and the ACT polarization 
(Naess et al. 2014). 

In order to describe the primordial tensor perturbations, we 
usually use the amplitude parameter, the tensor-to-scalar ratio, 
r = P T /P s , where P T and P s denote the power spectra of tensor 
and scalar perturbations. Due to the precision limitations of 
current CMB experiments, the constraints on r are still weak 
and consistent with zero very well (Montory et al. 2006; 
Hinshaw et al. 2013; Araujo et al. 2012; Kaufman et al. 
2013; Ade et al. 2013a, b). In 2014, the BICEP2 Collabora¬ 
tion published their final three-year CMB observational data 
and announced the detection of CMB B-modes polarization 
at scales 20 < £ < 340 (Ade et al. 2014), which corresponds 
to r ~ 0.2, and the null detection (r = 0) is disfavored at about 
la. 

However, few months later, the Planck Collaboration mea¬ 
sured the angular power spectra of the polarized thermal emis¬ 
sion from diffuse Galactic dust, which is the main foreground 
present in measurements of the CMB polarization at frequen¬ 
cies above 100 GHz (Gold et al. 201 1), and pointed out that 
the contribution of thermal dust on the CMB BB power spec¬ 
trum has been significantly underestimated in the BICEP2 
foreground analyses (Adam et al. 2014). They found that the 
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magnitude of the BB power spectrum, contributed by the po¬ 
larized thermal dust emission, at 150 GHz in the sky field ob¬ 
served by the BICEP2 experiment is at the same level as the 
one predicted by the early universe model with the primordial 
gravitational waves r ~ 0.2. When the polarized dust contri¬ 
bution is properly taken into account, the joint analysis of the 
BICEP2/Keck and Planck Collaborations preferred the null 
detection of the primordial gravitational waves and obtained 
the upper limit of the tensor-to-scalar ratio r < 0.11 at 95% 
confidence level (Ade et al. 2015a). 

Besides the polarized thermal dust emission, the cosmic 
birefringence can also generate the CMB B-modes polar¬ 
ization at large scales. This phenomenon can be introduced 
by the interaction between CMB photons and the external 
scalar field p t , through the Chern-Simons (CS) term £cs ~ 
p ll A„Fi- u/ . Here F pv = (1 /2)e^ pa F pa is the dual of the elec¬ 
tromagnetic tensor. This coupling is gauge invariant if d„p p = 
d p p v . This is possible if p p is a constant field over the 
spacetime or arises from the derivative of a cosmic scalar 
field (j). The scalar field 0 can be identified as the dark en¬ 
ergy in the quintessential baryo-/leptogenesis (Li et al. 2002; 
Li & Zhang 2003) and as the Ricci scalar R in the gravita¬ 
tional baryo-/leptogenesis (Li et al. 2004; Davoudiasl et al. 
2005). The CS term violates the Lorentz and CPT symmetries 
spontaneously if p f , has nonzero background values. This in¬ 
duces rotations of polarization directions of the propagating 
photons (Carroll et al. 1990; Carroll 1998) and productions 
of CMB TB and EB cross-correlations, even though they are 
absent before the recombination in the traditional CMB the¬ 
ory (Lue et al. 1999; Peng et al. 2005, 2006; Li et al. 2007; 
Xia et al. 2008a). The crucial point for this paper is that this 
phenomenon provides a mechanism to produce the CMB B- 
modes, alternative to the primordial gravitational waves and 
weak lensing. Because the rotation converts part of E-modes 
polarization to B-modes polarization, even the primordial B- 
modes polarization is absent, sizable CMB B-modes power 
spectrum can still be obtained from the E-modes power spec¬ 
trum through the rotation. As shown in Xia et al. (2010), the 
cosmological birefringence degenerates with primordial ten¬ 
sor mode perturbations at large scales and with the gravita¬ 
tional lensing at small scales, respectively. And due to this 
degeneracy, considering the anisotropic polarization rotation 
could slightly lower the best-fit value of r and relax the ten¬ 
sion on the constraints of r between from BICEP2 and from 
Planck 2013 data (Li et al. 2015), even without the help of 
the isotropic part of the rotation angle, which is removed 
by the “self-calibration” method in the BICEP2 data analy¬ 
sis (Ade et al. 2014). 

Therefore, it is important and necessary to study the influ¬ 
ences of the polarization rotation angle on the measurements 
of CMB primordial tensor perturbations using the current and 
future CMB experiments. Recently, Xia (2012) found that 
constraints on the tensor-to-scalar ratio could be significantly 
biased, when the isotropic rotation angle can not be taken into 
account properly. In this paper, our main task is to focus on 
the influences of the anisotropic polarization rotation angle on 
the detection of primordial gravitational waves by performing 
the global fitting with the latest CMB datasets and simulated 
mock data of future CMB experiments included. This paper 
is organized as follows: in section 2 we briefly introduce the 
cosmic birefringence mechanism and the effects on CMB po¬ 
larization power spectra from the isotropic and anisotropic po¬ 
larization rotation angle. Section 3 contains our main results 


from the current observations and future measurements, while 
section 4 is dedicated to the conclusions and discussions. 

2. CMB POLARIZATION ROTATION ANGLE 

In this section we briefly introduce the cosmic bire¬ 
fringence which has been widely studied in the litera¬ 
tures (Carroll et al. 1990; Carroll 1998; Lue et al. 1999; 
Peng et al. 2005, 2006; Li et al. 2007; Xia et al. 2008a; 
Li & Zhang 2008; Kamionkowski 2008; Yadav et al. 2009; 
Caldwell, Gluscevic & Kamionkowski 2011). We consider 
the CS coupling between the cosmic scalar field (j) and the 
electromagnetic field: 

C c& = d P M)Aj p \ (1) 

where / is an arbitrary function of 0. As mentioned be¬ 
fore, the CS term violates Lorentz and CPT symmetries but 
not the gauge symmetry of the electromagnetic field, so will 
not change the number of dynamical components of photons. 
However, the states with opposite helicity propagates with dif¬ 
ferent velocities, so that the polarization vectors of photons 
will not be parallel transported along the light rays. Conse¬ 
quently the observed polarization directions are changed by 
rotations when compared with the original directions at the 
source. 

In terms of the Stokes parameters Q and U, the rotated lin¬ 
ear polarization state (denoted by a prime) is related to the 
unrotated one through 

( Q' ± iU') = exp(±/2a)(e ± iU ). (2) 

The rotation angle a corresponds to the CPT violation and is 
frequency independent. It can be calculated by the integration 
along the light rays (Li & Zhang 2008) 

a= [ d p mdx p ( A) = /(&)-/(0o), O) 
J O 

where A is the affine parameter of the light ray, the subscripts 
s and o represent the source and the observer, respectively. 
Naturally, we observe CMB photons at a single point, hence 
c f> 0 should be a constant. However the dynamical field 0 is 
expected to fluctuate on the last scattering surface (LSS). The 
observer at the fixed position receives CMB photons from all 
directions and will find that the rotation angle varies across the 
sky, as pointed out in Li & Zhang (2008). The anisotropies 
are proportional to the distribution of <f> on the LSS, which 
means we can separate the rotation angle a into the isotropic 
part and the small anisotropic one as follows 

ct(n) = a + Aa(n), (4) 

with 

a = /[</>(r?LSs)]-/(</’o) , 

Aa(n)=^^(xLss,?? LSs ), (5) 

where ?/[ $s and j?lss denote the position and the conformal 
time on the last scattering surface. Lor the isotropic part 
a, the first evidence in terms of the full CMB datasets 
was done in Peng et al. (2006) and stimulated many inter¬ 
ests in this field (see Li et al. (2007); Xia et al. (2008a); 
Komatsu et al. (2009); Xia et al. (2008b); Wu et al. (2009); 
Brown et al. (2009); Komatsu et al. (2011); Xia et al. 
(2010); Liu et al. (2006); Xia et al. (2008c); Xia (2012); 
Hinshaw et al. (2013); Geng et al. (2007); Cabella et al. 
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(2007); Kostelecky & Mewes (2007); Kahniashvili et al. 
(2008); Finelli & Galaverni (2009); Li et al. (2009); 
Gruppuso et al. (2013); Zhao & Li (2013); Zhao et al. 
(2015), and references within). 

Assuming the anisotropic rotation angle Aa(n) is a zero 
mean field defined on two dimensional sphere, it can be de¬ 
composed in terms of spherical harmonics: 


Aa(n) = Y bimYtmi n) . (6) 

Furthermore, if we assume the perturbation field Aa(n) is a 
Gaussian random field and satisfies the statistical isotropy, so 
we can define the angular power spectrum of the anisotropic 
rotation angle: 

(b im b* ilm ,)=crsws mm ' ■ (7) 

In this paper we set f(cj)) = ap/M for simplicity, in which 
M represents the cut-off mass in the effective field theory 
and c stands for the dimensionless coupling coefficient. We 
find the relationship between the angular power spectrum of 
anisotropic rotation angle Cf a and the power spectrum of the 
scalar field P,p at the LSS 


<W> 


LSS 


2t r 2 
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P<f>(k)jj(kr]o-kr) s ) (9) 


where <j>^ is the Fourier transform of r)0 and ji(x) is the spher¬ 
ical Bessel function. 


J 


Zhao & Li (2014) has considered several specific dark 
energy models, such as the massless scalar field and the 
quintessence field, and then calculated their power spectra ex¬ 
plicitly. In this paper we study the constraint on the pertur¬ 
bations of this scalar field through the polarization rotation 
angle in a model-independent way. It is convenient to expand 
its power spectrum on the super hubble scale, just as what we 
usually do in the study of slow-roll inflation, 

/ k\ nceT 

P(j> = PCYY = ^CPT f ■£“ J , (10) 

where Acpt , «cpt and k {) represent the amplitude and spectral 
index of the power spectrum and the pivot scale, respectively. 
Due to the precision limitation of current CMB experiments, 
we set «cpt to be zero and neglect the higher order running 
parameters of the spectral index [dInncpr/d Ink] in the fol¬ 
lowing calculations. 

In order to investigate the influence on the CMB polariza¬ 
tion power spectra due to the cosmic polarization rotation an¬ 
gle, we firstly define two-point correlation function using the 
following relation 

2P +1 

C”03) = (fo(n)fo(n')} = Y Z ^CrPe( cos/3), (11) 

t 

where /3 corresponds to the angle between these two direc¬ 
tions, cos(/3) = n • n'. Combining Eqs. (2) and (6), one can ex¬ 
press the rotated CMB polarization power spectra C' e in terms 
of the unrotated ones Q (see Li & Yu (2013) for details): 


C ,EE +C' BB = exp [-4C°“(0)] Y (C EE + Cf) [ 4(/3)4(/3) e 4c ““ ( ^£/cos(^) 

cf E -C’ BB = cos(4a)exp[-4C a “(0)]^ -Cf) f d^ifod^py^^dc os(/3) 

C ,EB = sin(4a)exp[-4C“ Q (0)]y' ^-^(C^ £ -C^) [ 4 2 (/3)4-(/3)e“ 4C ““ (/3) ^cos(/3) 

_ o/'-f-l /*! t aa 

C'/ E = cos(2a) exp [~2C aa (0)] —-— Cj, E / d{ n ((3)d2 0 (f3)dcos((3) = C EE cos(2a)e~ lc (0) 

V ^ _I 

2 p f -L. i /*i t 

C ,EB = sin(2a)exp [~2C aa (0)] ^ A_— cJ, E / dQ 2 (P)d 2 o(P)dcos(p) = C EE sin(2a)e~ 2c (0) , 

«/ ^ 


( 12 ) 


where C a “(0) = ’Yhp(2(.+ l)C““/47r is the auto-correlation of 
the anisotropic polarization rotation angle. Note that, in the 
calculations we have assumed the original CMB fields and 
anisotropic rotation angle field are statistically isotropic Gaus¬ 
sian random fields. 

In the Eqs. (12) we have separated the influences of 
isotropic and anisotropic polarization rotation angle. It is ap¬ 
parently that they both get involved in all scales of CMB po¬ 
larization power spectra. The isotropic rotation angle copies 
the original CMB EE power spectrum, zooms it out and 
then prints it onto the rotated BB power spectrum. It will 
inevitably create additional CMB TB and EB power spec¬ 
tra which should be vanished in normal CMB theories. We 


can see the terms cos(a) or sin(a) appear as coefficients in 
almost all equations. This suggests us that the isotropic ro¬ 
tation angle should perform a global modulation on rotated 
CMB power spectra if it exists. However, the situation for 
the anisotropic one is somehow different. Non-zero TB and 
EB spectra are not longer necessary if the isotropic part [a] 
vanishes, which is consistent with some current CMB ex¬ 
periments. The BICEP2/Keck Array collaboration has used 
the “self-calibration” method for the detector polarization ori¬ 
entations, and as a result, any signal of the isotropic polar¬ 
ization rotation angle has been eliminated from data. Fortu¬ 
nately, the anisotropic part is not affected by this calibration 
method, and will still contaminate the observations of the pri- 
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Fig. 1.— The theoretical CMB BB power spectra for different values of 
1 n[ 1 0 4 .4 cptI- We also show the theoretical CMB BB power spectrum with 
the tensor-to-scalar ratio r = 0.02 (magenta thin line) for comparison. 


mordial gravitational waves. In figure 1, we plot the CMB 
BB power spectra for different values of 1ii[10 4 Acpt] with 
a = 0. For comparison, we also show the expected CMB BB 
power spectrum generated by the primordial tensor perturba¬ 
tions r = 0.02. Interestingly, we find that, even the isotropic 
rotation angle is zero, the anisotropic polarization rotation an¬ 
gle with ln[10 4 Acpr] ~ 2.5 could still generate the CMB BB 
power spectrum to mimic CMB primordial B-modes signal 
with r = 0.02 at large scales. At scales £ < 20, the contribution 
on BB power spectrum from the anisotropic polarization ro¬ 
tational angle can be negligible. Therefore, probing the reion¬ 
ization bump of CMB BB power spectrum at very large scale 
may provide a smoking gun for distinguishing the polariza¬ 
tion rotational angle and primordial gravitational waves. 

3. NUMERICAL RESULTS 

In this paper we perform a global analysis to all the power 
spectra of the CMB data with the public available Markov 
Chain Monte Carlo package CosmoMC (Lewis & Bridle 
2002), which has been modified to compute the rotated CMB 
polarization power spectra discussed above. We assume the 
purely adiabatic initial conditions and impose the flatness con¬ 
dition motivated by Inflation. Our basic parameter space is: 
P = (< jJb,u) c ,Q.\,T,n s ,A s ,r ), where u>b = Clbh 2 and uj c = Q. c h 2 
are the physical baryon and cold dark matter densities relative 
to the critical density, U \ is the dark energy density relative to 
the critical density, t is the optical depth to re-ionization, A s 
and n s characterize the primordial scalar power spectrum, r is 
the tensor to scalar ratio of the primordial spectrum. For the 
pivot of the primordial spectrum we set k s o = 0.05 Mpc -1 . For 
the polarization rotation angle, we set the isotropic rotation 
angle a = 0 and then have one free parameter, Acpt, which 
represent the amplitude of power spectrum of the anisotropic 
rotation angle. Furthermore, in our analysis we include the 
CMB lensing effect, which also produces B-modes from E- 
modes (Zaldarriaga & Seljak 1998), when we calculate the 
theoretical CMB power spectra. 

3.1. Current CMB Observations 

In our calculations we mainly use the polarization power 
spectra from BICEP2 experiments and the joint analysis data 
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Fig. 2.— One- and two-dimensional constraints on the tensor-to-scalar ratio 
and the anisotropic polarization rotation angle from the BICEP2 data (black 
solid lines). For comparison, we also show the constraint on r from the BI- 
CEP2 data with Acpt = 0 (red dashed line). 


of BICEP2/Keck Array and Planck experiments: 

• The Background Imaging of Cosmic Extragalactic Po¬ 
larization (BICEP2) (Ade et al. 2014), locating at the 
South Pole, specializes in searching polarized CMB 
signal at large scales with low angular resolution, 
concentrating 150 GHz on a roughly 1% patch of 
sky at high Galactic latitudes. Recently the BICEP2 
collaboration released their three-year data accumu¬ 
lated from 2010 to 2012 and announced the detec¬ 
tion of CMB B-modes polarization at scales 20 < 
£ < 340. It is noteworthy that they claimed that the 
CMB EB power spectrum is only used for the “self¬ 
calibration” of the detector polarization orientations 
(Keating, Shimon & Yadav 2013). Any polarization 
rotation has been removed from the results. 

• The products of the Joint Analysis of BICEP2/Keck Ar¬ 
ray and Planck Data (BKP2015) (Ade et al. 2015a). 
Very recently, BICEP2/Keck Array and Planck Col¬ 
laborations have re-analysed their data and published 
the CMB BB power spectrum at scales 20 < £ < 
340, in order to investigate the influence on the CMB 
B-modes measurements from the polarized thermal 
dust emission. They convert the Planck maps into the 
format usable for BICEP2/Keck Array experiments, 
and make cross-spectra between the combined BI- 
CEP2/Keck maps and Planck maps of the all polar¬ 
ized bands. Similar with the BICEP2 data, the “self¬ 
calibration” method was also applied in the data analy¬ 
sis. 

Due to the “self-calibration” method used in the BI- 
CEP2/Keck Array experiments, the constraints on the 
isotropic rotation angle [a] are quite consistent with zero 
(Li et al. 2015). Therefore, we only use these data to study 
the anisotropies of the polarization rotation angle. Further¬ 
more, in order to take the contribution of polarized thermal 
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Fig. 3.— One-dimensional constraints on the tensor-to-scalar ratio from 
the BKP2015 data with (black solid line) and without (red dashed line) the 
effect of the anisotropic polarization rotation angle included. 



r 


Fig. 4.— One-dimensional constraints on the tensor-to-scalar ratio from 
the simulated mock CMB data with (black solid line) and without (red dashed 
line) the effect of the anisotropic polarization rotation angle included. 


dust emission on the CMB B-modes power spectrum into ac¬ 
count properly, we add one more nuisance parameters ,4 ( i us i 
in our calculations and consider the polarized dust BB power 
spectrum satisfies C® dust oc Ad ust £- 2A2 (Adam et al. 2014). 

In figure 2 we show the constraints on the tensor-to-scalar 
ratio and the anisotropic polarization rotation angle, and their 
2-dimensional contours from the BICEP2 data. Firstly, we ne¬ 
glect the contribution on the CMB B-modes power spectrum 
from the polarization rotation angle, which means Acpt = 0. 
Due to the strong degeneracy between r and ,4 dlls t, the huge 
uncertainty of the polarized dust B-modes power spectrum 
significantly reduces the constraining power of BICEP2 data 
on the tensor-to-scalar ratio r. The obtained 2cr upper limit is 
r < 0.12 and there is a peak around r ~ 0.05, shown as the red 
dashed line in figure 2. There is no signal for the primordial 
gravitational waves, which is consistent with the results from 
the BKP2015 data (Ade et al. 2015a). 

When we include the effect of anisotropic polarization rota¬ 
tion angle in the calculations, the constraint on the tensor-to- 
scalar ratio is slightly different. As we mentioned in Li et al. 
(2015), the anisotropic polarization rotation angle would also 
contribute to the CMB BB power spectrum and partly ex¬ 
plain the CMB B-modes data of BICEP2. In this case, the 
1-dimensional distribution of r slightly moves towards to the 
smaller value, due to this degeneracy, see the black solid line 
in figure 2. And the minimal value of x 2 becomes slightly 
lower. Ax 2 ~ — 1. We also check this result using the lat¬ 
est BKP2015 data, which is shown in figure 3. When we 
include the contribution on CMB B-modes power spectrum 
from the anisotropic polarization rotation angle, the upper 
limit of tensor-to-scalar ratio is r < 0.13 (95% C.L.) and 
the peak structure totally disappears. These results imply that 
even considering the effect of polarized thermal dust emis¬ 
sion, the anisotropic polarization rotation angle can still affect 
the constraint on the tensor-to-scalar ratio slightly. 

In figure 2 we also show the one-dimensional constraint 
on the amplitude parameter Acpt- The current BICEP2 data 
still do not have enough constraining power to limit Acpt 
accurately, so the obtained constraint on Acpt is consistent 
with zero, namely the 95% upper limit is ln[10 4 AcFr] < 3.8. 
However, we interestingly find that there is a peak around 


1ii[10 4 Acpt] ~ 2.5. As we show in figure 1, the model with 
ln[10 4 AcpT] ~ 2.5 could mimic the contribution of CMB B- 
modes power spectrum of primordial gravitational waves with 
r = 0.02 at large scales. And this degeneracy moves the distri¬ 
bution peak of r towards to the lower value, which is similar 
with the results in Li et al. (2015). 

Finally, we plot the two-dimensional contour on the panel 
(r,ln| 1 0 4 4cpi]) in figure 2. Apparently, these two parame¬ 
ters are slightly correlated. The model with a larger value of 
ln[10 4 AcpT l provides more power of CMB B-modes, and con¬ 
sequently the contribution from the primordial gravitational 
waves is suppressed, namely the model with r = 0 is more 
favored by the BICEP2 data. Conversely, if we limit the pa¬ 
rameter ln[10 4 Acpr] to be a very small value, the large value 
of r will be needed, in order to match the BICEP2 data. 

3.2. Future CMB Measurements 

Since the current CMB data are not accurate enough, the 
degeneracy between r and ln[10 4 AcpT] is not very clearly. 
But this could be important in the future CMB B-modes 
experiments. Therefore, it is worthwhile discussing the ef¬ 
fect of non-zero anistropic polarization rotation angle on the 
measurement of CMB primordial B-modes signal for the 
future experiments. Following Xia et al. (2008c), we sim¬ 
ulate the future mock CMB power spectra with BICEP3- 
like sensitivity: /Ay = 0.05 sky converge, the isotropic noise 
with variance At = Ap/2 = 0.13/iX 2 , a symmetric Gaussian 
beam of 25 arcminutes full-width half-maximum (FWHM) 
(Ahmed et al. 2014), and the maximum multipoles £ max = 
1000. Here, we neglect the systematic errors of future CMB 
measurements and the contaminations from the CMB vari¬ 
ous foregrounds, such as the thermal dust emission and the 
synchrotron emission. The basic fiducial model we choose is 
the best-fit Planck model (Ade et al. 2015b): tti,h 2 = 0.02225, 
Q c h 2 = 0.1194, 1000, = 1.04094, r = 0.079, n s = 0.9682, 
log [10 10 A S ] = 3.08919 at k s0 = 0.002Mpc' 1 . 

Firstly, we simulate the mock data using the basic fidu¬ 
cial model and the tensor-to-scalar ratio r = 0.03. When us¬ 
ing the basic seven cosmological parameters (Acpt = 0) to fit 
the mock data, we obtain a 3cr detection of non-zero value 
of tensor-to-scalar ratio: r = 0.03 ±0.01 (68% C.L.). How- 
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Fig. 5.— The one-dimensional posterior distributions of the tensor-to- 
scalar ratio r derived from the simulated mock CMB data without includ¬ 
ing the anisotropic rotation angle. The anisotropic rotation angles consid¬ 
ered in the fiducial model are: ln[10 4 Acpj] = 0 (blue dash-dotted line), 
ln[10 4 AcpTl = 1 (red dashed line), and ln[10 4 AcpT] = 2 (black solid line). 


ever, when we set the parameter ln[10 4 AcpT] to be free, the 
constraint on r is obviously weakened, r = 0.03 ± 0.03 (68% 
C.L.) and r < 0.1 (95% C.L.), as shown in figure 4. Note that, 
the standard deviations of r we obtain are apparently under¬ 
estimated, since in the calculations we do not consider any 
contamination from the CMB foregrounds. This result implies 
that in the future CMB experiments, if we want to measure 
the signal of CMB primordial B-modes accurately, we need 
to know very well about the anisotropic polarization rotation 
angle. Otherwise, the degeneracy among them will contami¬ 
nate the signal of primordial gravitational waves significantly. 

We also simulate another mock data using the basic fidu¬ 
cial model and the non-zero anisotropic polarization rotation 
angle. Here, we consider three fiducial values ln[10 4 AcFr] = 
0,1,2. When we include the parameter of anisotropic polar¬ 
ization rotation angle in the calculations, the fiducial values 
of input parameters are always recovered. However, if we ne¬ 
glect the effect of anisotropic rotation angle at all Acpt = 0, 
the contribution of non-zero anisotropic rotation angle on the 
CMB BB power spectrum will be wrongly considered as the 
signal of the primordial CMB tensor perturbations. Since the 
fiducial model r - 0, the BB power spectrum due to the pri¬ 
mordial tensor B-modes should vanish. But in the presence 
of non-zero fiducial rotation angle (ln[10 4 AcpT] = 0,1,2), the 
BB power spectrum should be non-vanishing. If we force the 
anistropic rotation angle to be zero in the analysis, the value 
of r will be enlarged to match the mock non-zero BB power 
spectrum. In figure 5, we find that the larger the fiducial value 
ofln[10 4 A CPT ] is, the larger the obtained central value of r be¬ 
comes. For example, in the calculation with the fiducial model 
ln[10 4 AcpTl = 2, if we fix Acpt = 0, the obtained constraint on 
the primordial tensor perturbations r = 0.017 ±0.006 (68% 
C.L.). The significance of this fake signal is about 3cr confi¬ 
dence level. Therefore, if we do not take the effect of polar¬ 
ization rotation angle into account properly, the constraints on 
the tensor-to-scalar ratio r from the future CMB data with the 
high precision will be significantly biased. 

4. CONCLUSIONS AND DISCUSSIONS 


One of the most important scientific aims in the CMB cos¬ 
mology is searching for the primordial gravitational waves on 
the CMB B-modes power spectrum. However, the signal is 
very small and could be contaminated by other sources easily, 
such as the CMB foregrounds, polarized thermal dust emis¬ 
sion and synchrotron emission. In this paper, we discuss an¬ 
other mechanism, the cosmic birefringence, which can be in¬ 
troduced by the CPT-violating interaction between CMB pho¬ 
tons and the external scalar field. This mechanism can also 
produce the CMB B-modes power spectrum and mimic the 
signal of primordial gravitational waves at large scales. Here 
we summarize our main conclusions in more detail: 

• We use the BICEP2 data and the joint analysis data of 
BICEP2/Keck Array and Planck experiments to con¬ 
strain the tensor-to-scalar ratio. When we include the 
contribution of polarized thermal dust emission on the 
CMB B-modes power spectrum, the constraints on r are 
consistent with zero at 95% confidence level. Further¬ 
more, when we consider the effect of the anisotropic 
polarization rotation angle into the analyses, the degen¬ 
eracy between r and Acpt slightly changes the con¬ 
straint on r. The one-dimensional distribution of r 
slightly moves towards to the smaller value, since the 
anisotropic polarization rotation angle would also con¬ 
tribute to the CMB BB power spectrum and partly ex¬ 
plain the CMB B-modes data. 

• The current BICEP2 data still do not have enough con¬ 
straining power to study the polarization anisotropic ro¬ 
tation angle accurately, so the obtained constraint on 
Acpt is consistent with zero. However, there is a clear 
peak around ln[10 4 AcpT] ~ 2.5, which could mimic the 
contribution of CMB B-modes power spectrum of pri¬ 
mordial gravitational waves at large scales, as shown in 
figure 1. We also find the slight correlation between r 
and ln[10 4 AcpT] in the two-dimensional contours. 

• We simulate the mock CMB data using the BICEP3- 
like sensitivity and study the degeneracy between the 
tensor-to-scalar ratio and the anisotropic rotation angle. 
When we neglect the effect of nonzero Acpt, the mock 
data can constrain r tightly. However, when including 
this effect into analysis, the obtained standard deviation 
of r is significantly enlarged by a factor of three. This 
degeneracy will contaminate the signal of primordial 
gravitational waves obviously. 

• Finally, we simulate the mock data with the different 
fiducial values. If we neglect the effect of anisotropic 
rotation angle in the calculations, the contribution of 
non-zero anisotropic rotation angle on the CMB BB 
power spectrum will be wrongly considered as the sig¬ 
nal of the primordial tensor perturbations. Therefore, 
we need to take the effect of polarization rotation angle 
into account properly. Otherwise, the constraints on the 
tensor-to-scalar ratio r could be biased from the future 
CMB experiments with high precisions. 
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